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ABSTRACT 

We present a weak lensing analysis of the Coma Cluster using the Sloan Dig- 
ital Sky Survey (SDSS) Data Release Five. Complete imaging of a ~ 200 square 
degree region is used to measure the tangential shear of this cluster. The shear is 
fit to an NFW model and we find a virial radius of r2oo = l-99ln22h -1 Mpc which 



corresponds to a virial mass of M 2 oo = 1.88io!56 x 10 15 h _1 M Q . We additionally 
compare our weak lensing measurement to the virial mass derived using dynami- 
cal techniques, and find they are in agreement. This is the lowest redshift, largest 
angle weak lensing measurement of an individual cluster to date. 
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Subject headings: galaxies: clusters: individual (Coma) — gravitational lensing 



1. Introduction 



Measurements of the mass of the Coma cluster date back to IZwickyl (119331 ) who first 
applied the virial theorem to Coma and showed that dark matter dominates the cluster on 
megaparsec scales. Since then a num ber of methods have been used to determine the mass 
of Coma under various assumptions. Kent fc Gunnl (119821 ) used the galaxy surface density 
profile t o measure the mass under the assumption that the mass traces the galaxy distri- 
bution. The fc White! (119861 ) l ater use d this dataset to measure the mass with a modified 
version of the virial theorem. iHughesI (119891 ) used X-r ay data und e r the assumption that 
the cluster is in hydrostatic equilibrium. More recently iGeller et al.l (119991 ) used the cluster 
infall region of Coma to extrapolate the mass out to 10/i -1 Mpc as well as to determine the 
virial mass. 
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Wea k lensing has beco me a powerful tool with which to measure the mass of galaxy 
clusters (ITyson et al.lll990l ). Here the shape distortion of background galaxies due to a 
foreground cluster is used to determine the mass of the cluster independent of dynamical 
assumptions. The shear induced on the background galaxies can be used to probe the 
mass of the foreground cluster out to large radii. The m ajority of galaxy clusters studied 



with weak lensing lie at redshift z > 0.2 (IDahk- 120071). However lower redshift clusters 



have been probed using this method, for instance iJoffre et al.l ( 120001 ) studied Abell 3667 at 
z = 0.055. At very low redshifts weak lensing measurements of clusters become difficult 
since this requires imaging of a large area surrounding the cluster. 

The possibility of a weak lensing study of the Com a Cluster (z = .0236 ) with the SDSS 
was first suggested by I Gould Villumsenl (11994 ) and lStebbins et al. Jl996h . They pointed 
out that since Coma is at such a low redshift the weak lensing signal would be measurable 
in the SDSS since virtually all background galaxies would be at a much higher redshift than 
Coma. This is in spite of the shallow imaging in the SDSS which yiel ds a typical background 
galaxy surface density of ~ 1 galaxy arcmin -2 ( jSheldon et al.l 120041 ) . Complete imaging of 
the entire reg ion surrounding the Coma Clus ter was recently completed in the SDSS Data 
Release Five (jAdelman-McCarthv et al.ll2007l ). 



Here we present a weak lensing measurement of the virial mass of the Coma Cluster 
in the SDSS. Our paper is organized as follows: in §2] we discuss our weak lensing analysis, 
in §3] our determination of the mass of Coma is discussed, and in £JH we compare our weak 
lensing mass estimate to previous measurements of the virial mass of Coma. 

Since Coma is at such a low redshift our results do not change significantly for different 
density parameters. However for completeness we have assumed a standard cosmology with 



0.3 and ft A = 0.7. 



2. Weak Lensing Analysis 



2.1. Data 



Data used in our study are obtained from the Sloan Digital Sky Survey (SDSS) (jYork et al. 



20001 ) a large imaging and spectroscopic survey of an 8000 square degree region in the 
Northern Galactic Cap centered on a = 12 h 22 m , S = 32° 13' (J2000). The SDSS uses a 
dedicat ed 2.5m telescope at Apache Point Observ atory which images th e sky in the ugriz 



bands (IFukugita et al.lll996l ) in a drift scan mode 



Tucker et al. 


(2006 


) and 


Hogg et al. 


(2001) 



Stoughton et al.l 120021 ). The astrometric 
id the photometric calibration is 



- 3- 



Data Release Five (lAdelman-McCarthy et al.ll2007l ) which contains complete imaging of the 
region surrounding Coma. 

Object detectio n and shape measurement are performed using the PHOTO pipeline 
( iLupton et al.l 12001). PHOTO measures a large number of parameter s for each detected 
objec t in the SDSS (jStoughton et al.ll2002l ). including adaptive moments ([Bernstein fc Jarvis 
20021 ) which we use to measure galaxy shapes. In the SDSS adaptive moments are measured 
using an iter ative algorithm . which adapts a Gaussian weight function to the size and shape of 
each galaxy (ISheldon et al.ll2004l ). PHOTO convolves each imag e with the local PSF befor e 
detecting objects in order to avoid the selection bias described in lBernstein fc Jarvid (120021 ). 
Drift scanning in the SDSS creates a time dependent PSF which leads to a spatial variation 
in image quality. To model the PS F in the SPSS a Kar hune n-Loeve (KL) d ecomp osition is 
used, described in further detail in iLupton et al.l (120011 ) and ISheldon et al.l (120041) . Camera 



shear due to drift scanning in the SDSS has been previously shown to be small (IHirata et al. 
20041 ). 



For our weak lensing analysis we use a ~ 200 square degree region centered on the 
core of the Coma Cluster. The core of Coma contains two bright cluster galaxies (BCGs) : 
NGC 4874 & NGC 4889. We take as the center of Coma the galaxy NCG 4874 located at 
a = 13 h 02 m 00.2 s , 5 = 2 7°4l'26.6" (72000). T his BCG is associated with the majority of the 
X-ray flux in the core (jVikhlinin et al.ll200ll ) and the large number of galaxies surrounding 
this BC G in the optical imag ing is indicative of this being the center of the cluster potential 
(see e.g., lHansen et al.ll2005l ). We in fact find that the shear amplitude is maximized when 



we use this BCG as the cluster center. 



2.2. Source Galaxies 



To create our source galaxy catalog we use only objects which have been detected in all 
three gri filters. We eliminate objects which contain saturated pixels, were originally blended 
with another object, or triggered errors in the measurement of adaptive moments by rejecting 
objects with the respective PHOTO error flag set. For our study we use only objects which 
have been classified by PHOTO as galaxies (type = 3). We additionally use only shape 
measurements from the r band since this filter is the most sensitive and typically contains 
better seeing jAdelman-McCarthv eTall Ig^oj)^ Galaxies are selected in the r band with 
extinction corrected model magnitudes (jStoughton et al.l 120021 ) in the range 18 < r < 21 
(Figured]). Model ma gnitudes are used here in stead of Petrosian magnitudes since these 
provide a higher S/N (IMandelbaum et al.l 120051 ). Galaxies with fainter magnit udes have 
been used in galaxy-galaxy weak lensing studies in the SDSS (ISheldon et al.ll2004l ). but here 
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we restrict our sample to the ma gnitude range where errors in shape measurements are 
typically small (IHirata et al.ll2004l ) and galaxy sizes are larger. 



In our analysis we us e photometric redshifts for our source galaxies. The SDSS photo- 
metric redshift pipeline of ICsabai et al.l (120031 ) is used since it completely covers the area of 
Coma in Data Release Five. Galaxies in our magnitude range are typically at much higher 
redshift than Coma, but we use photometric redshifts here to aid in the rejection of fainter 
cluster members. For our study we use only galaxies with photometric redshifts in the range 
0.2 < Zphot < 0.8 and redshift errors z err < 0.4. Galaxies with lower photometric redshifts 
are not used as the fractional error in redshift rapidly increases and most of the shear signal 
comes from galaxies at high redshifts compared to Coma. 

To correct for t he PS F anisotropy and dilution we use the linear PSF correction scheme 
of IHirata fc Seljakl (120031 ). This algorithm uses the measured galaxy ellipticity and recon- 
structed PSF at the position of the object to correct the ellipticity components of each galaxy. 
We use galaxies with corrected elli pticities e mr r < 1. 4 in our analysis, which is typical for 
weak lensing studies in the SDSS (IHirata et al.ll2004l ). We additionally make a cut on the 
resolution parameter (R), defined as 



R 



M PSF 

j rrcc 

fvl rrcc 



(IBernstein fc Jarvid 120021 ) where M rrcc is the size of the galaxy and is the size of the 

PSF at the position of the galaxy as calculated by PHOTO. For our analysis we use only 
galaxies with R > 0.33, which is equivalent to only using objects 1.5 times larger than the 
PSF. We additionally rotate the PSF corrected ellipticity components of each galaxy from 
image coordinates to the equatorial coordinate system. 

Within the annulus described in §2.31 the total number of galaxies in our source galaxy 
catalog after applying these cuts is ~ 270, 000. 



2.3. Shear Measurement 

The shear due to Coma is measured by projecting the PSF corrected ellipticity com- 
ponents of each source galaxy to the tangential frame and binning the galaxies into radial 
annuli between 0.05h^ 1 Mpc and 10.5h _1 Mpc. The edge of the inner annulus is chosen to 
be larger than the Einstein radius of Coma (~ 30") to ensure we are outside of the strong 
lensing regime, and also to avoid contamination from the BCGs. Since our radial annuli 
extend out to ~ 10h _1 Mpc (~ 8°), we modify the usual flat sky shear matrix into a curved 
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sky shear matrix, following ICastro et al. fl2005h . The tangential shear is measured using 

It = — ^ (2) 

where 74 is the tangential shear, et is the tangential ellipticity, N is t he total number of 
object s in each radial bin, and 1Z is the shear responsivity described in iBernstein fc Jarvis 
( 120021 ) . We chose not to use any ellipticity error weighting here since our galaxies are selected 
i n a magnitude ran ge where source galaxy ellipticity error is typically small in the SDSS 
( IHirata et al.ll2004l ). We additionally also chose not to weight each source galaxy by the 
critical surface mass density because all of our source galaxies are at a much higher redshift 
than Coma and therefore the weights are very nearly equal. In the case of no weights, the 
shear responsivity is 1Z — 1 — Og N , where <tsn is the shape noise (the width of the intrinsic 



ellipticity distribution per ellipticity compon ent') f|Wittmanll20 02'). Here we use a shape noise 
of cgN = 0.37 which was used previously in IHirata et al.l (120041 ) for the same source galaxy 
magnitude range. 

The resulting shear profile for Coma is shown in Figure [2j We expect to detect the shear 
out to ~ 5h _1 Mpc, and the shear (solid squares) is clearly detected in these inner bins. Also 
shown in Figure [2] is a null test (open triangles) where each source galaxy is rotated by 45°. 
The shear signal should dissappear after this rotation if the signal is due to lensing. Errors 
shown for both quantities in Figure [2] are la error bars, where a is the standard deviation 
of the mean in each radial bin. 

To statistically test the significance of our measured shear signal we test whether it 
is consistent with the null shear model (7$ = 0). We find that our measured shear gives 



a X 



23.33 for 6 degrees of freedom for the null model. The probability of obtaining a 



reduced \ 2 greater than this value is P(x 2 > 



xl) 



0.07%, and therefore we can reject the 



null hypothesis. For comparison, fitting the shear from the 45° test to the null model gives 
a x 2 — 5-65 for 6 degrees of freedom, which is consistent with the null hypothesis. 



2.4. Blank Fields Test 

To test that our signal is not affected by any remaining systematic in the survey, we 
split up the SDSS North into separate, non-overlapping 'Coma sized' patches that contained 
no galaxies used in our Coma analysis and which contained no large sections of missing data. 
We were able to successfully extract 6 blank fields in DR5 each ~ 20° x 20° wide. From the 
center of each patch we probed radially outward to ~ 10h _1 Mpc at the redshift of Coma. We 
applied the same cuts to the source galaxy catalog in each of the blank fields that were used 
in our Coma analysis in §2.21 The resulting inverse variance weighted average signal over all 
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fields is shown in Figure [31 where we have used the same binning as in our Coma analysis. 
The tangential shear signal is shown with solid squares and the 45° component is shown with 
open triangles. The tangential shear component is consistent with the null model giving a 
X 2 = 1.84 for 6 degrees of freedom. The 45° component gives a x 2 = 12.58 for 6 degrees of 
freedom which has a probability P(x 2 > Xo) = 5.0%. A model with a small positive shear 
is slightly favored for the 45° component, however the null model can only be excluded at 
the 5% level. The field to field scatter in the blank field shear measurements shows a typical 
standard deviation in a given radial bin of cx 7 ~ 0.0015 or a standard deviation of the mean 
of ~ 0.0006, similar in magnitude to the statistical errors that are plotted in Figure O In 
principle the sca tter over the bl ank fields can provide information on the error due to large 
scale structure (jHoekstral 1200 It also §4.41 below), but our determination of the scatter is 
limited by the finite number of blanks fields. 



Mass Model 



3.1. Tangential Shear 

The tangential shear due to a foreground cluster lens is given by 

£(< r) - S(r) 



It 



(3) 



crit 



where £(< r) is the average projec ted surface mass den sity interior to r, and S(r) is the 
projected surface mass density at r (IMiralda-Escudd 1 1 99 ll ) . The magnitude of the shear also 
depends on the critical surface mass density E crit which is determined by 



J crit 



AnG DiDi s 



(4) 



Here D\ and D s are the angular diameter distances from the observer to the lens and source 
respectively, and Di s is the angular diameter distance between the lens and source. 

To compute the critical surface mass density we use for the lens the exact spectro- 
scopic redshift of Coma z = 0.0236 (IGeller et al.lll999l ). We assume here that the peculiar 
velocity for Coma is zero, which has been shown in several studies ( IScodeggio et al.l 119971 ; 
Giovanelli et al.lll997l ). Photometric redshifts are used to calculate the source angular diame- 
ter distances, and we obtain a critical surface mass density of 26299M©pc~ 2 . For comparison 
this is nearly the value obtained if we assume all of our sources were at a fixed redshift of 
z = 0.3. 
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3.2. NFW Profile 



Since the S/N of our data is low (S/N ~ 5) we chose to fit the measured shear profile 
of Coma to a model. We fit the sh ear to that expected from a Navarro, Frenk, & White 
(NFW) profile (INavarro et al.lll996l ). The NFW profile is a "universal profile" found in N- 
body simulations to fit mass density profiles ranging from galaxies to galaxy clusters. The 
density of an NFW profile is described by 

S c Pc 



p(r) 



(5) 



(r/r a )(l + r/r s ) 2 

where 5 C is the halo overdensity, r s is the scale radius, and p c = 3H 2 (z)/8tiG is the critical 
density at the redshift of the cluster. The halo overdensity is given by 

200 I 
3 ln(l + c) - c/(l + c) 1 ' 

where c is the halo concentration. The NFW model is therefore determined by only two 
parameters c and r s which are highly correlated. The NFW model can be used to determine 
the virial radius r 2 oo = cr s defined as the radius in which the interior mass density falls to 
200p c . The corresponding virial mass M 20 o is given by 

8007T 



M- 



200 



3 

~Pc r 200- 



(7) 



The expected shear due to an NFW profile has been worked out in detail by lWright fc Brainerd 



( 120001 ) and we use their result in our analysis. Because the NFW profile is highly non-linear, 
we use the Levenberg-Marquardt fitting procedure in iPress et al.l (119951 ) to determine the 



virial radius and the halo concentration. We find that the values which minimize x 2 are 
given by 



r 200 



l^t^h^Mpc 



and 



3.84 



-13.16 
-1.84 • 



(9) 

The concentration parameter is not well constrained. This fit gives a x 2 — 3.87 for 4 degrees 
of freedom which has a P(x 2 — Xo) — 42.4% probability of occuring. This corresponds to a 
virial mass of 

M 200 = 1.88i°;| x 1O 15 M . (10) 

Our values of r 2 oo and M 2 oo remain constant if the number of bins is slightly changed, and 
therefore our choice of binning does not bias this result. It has been shown previously that 
l arge scale struc ture affects the errors in the fit parameters of weak lensing measurements 
(jHoekstral 120031 ) . however this is not considered in our fit errors. We discuss additional 
sources of systematic error not included in our fit in §4.41 
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4. Discussion 

We compare our weak lensing mass estimate with three other methods used to determine 
the virial mass of Coma: (1) The Virial Theorem (2) X-ray and (3) the Infall Region. Our 
value for M 2 oo is consistently higher than the previous measurements, though never by more 
that 2a. These comparisons are summarized in Table [TJ We also discuss sources of additional 
error not included in our mass measurement, as well as the source of any potential bias. 



4.1. Virial Theorem 



The fc White! (119861 ) used a sample of galaxies in Coma along with the virial theorem 
to measure its mass. Their analysis also accounted for an additional portion of the system 
that may have been missing in their spectroscopic sample. Relaxing the assumption that 
light traces mass in the system, they found that the mass varied over a wide range. However, 
assuming that the mass traces the galaxy distribution they found a mass within 2.7h" 1 Mpc 
of 0.95 x lO 15 h _1 M with a 15% error in mass. Our weak lensing measurement is consistent 
with this result to within ~ 1.6cr. 



4.2. X-ray 



Ou r weak lensing measurement is also consistent with the previous X-ray analysis of 



Hughesl (119891 ). This study used the assumption of hydrostatic equilibrium but avoided 



certain assumptions about the state of the gas and also accounted for many systematic 
effects. The best fit model was one in which mass traces light, giving a virial mass of 
(0.93 ± 0.12) x lO 15 h _1 M within 2.5h" 1 Mpc which is within - 1.7a of our result. 



4.3. Infall Region 



Using a redshift survey of ~ 1000 galaxies iGeller et al.l (119991 ) used the infall region 
of Coma to determine its mass. In redshift space galaxies falling into the potential well 
of the cluster are used to determine the amplit ude of the redshift caus tics (the boundaries 
in line of sight velocity vs. projected radius) (jDiaferio fc Gellerl 119971 ) . The amplitude of 
the caustics along with the assumptions of spherical symmetry and hierarchal clustering 
are directly related to the cluster gravitational potential. Using this technique they find 
an NFW model fits the mass well giving an estimate for r2oo = 1.5h -1 Mpc which yields an 
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^200 = 0.8 x 1O 15 M . Error bars are not quoted for their virial radius or mass; however, 
our weak lensing virial mass is in agreement with their central result to within ~ 2.0cr 



4.4. Additional Error and Bias 



Weak lensing measurements of the mass of clusters, in particular low redshift clusters, 
suffer from two sources of a dditional error: (1) distant large scale structure and (2) correlated 
structure near the cluster jHoekstralboolj ). Weak lensing is sensitive to all structure along 
the sight and it has been shown that any ba ckground struc ture introduces additional error 
in the weak len sing mass estima te of clusters (jHoekstrall200ll ) . Specifically for the case of the 
Coma Cluster, iHoekstral (120011 ) showed that the error due to large scale structure increases 
with imaging depth, thus limiting the total achievable S/N for Coma to ~ 7. In our study 
the imaging depth in SDSS is relatively shallow so the contribution of large scale structure 
to the total error budget is small. For future deep ima ging surveys which will also cover the 
area of Coma such as PanSTARRS (IKaiser et al.ll2002l ). the error due to large scale structure 
will be more significant. 

The effect of correlated structure (such as fi laments) ne ar a foreground clu s ter ha s been 
previously studied using N-body simulations by ICenl (119971 ) and Metzler et al.l (11999I ). Here 
it was shown that filaments can increase the statistical error in weak lensing mass estimates 
as well as cause the measured M 2 oo to be biased toward higher values. In principle Coma 
should be an interesting testing ground to study this effect since the neighboring structure is 
well understood. However in practice there is not a clear method to correct a weak lensing 
measured value of M 2 oo for this effect, and therefore we leave a study of any potential bias 
in our measurement for a future paper. 



5. Conclusion 



We have measured the weak lensing shear due to the Coma Cluster, currently the lowest 
redshift (z = 0.0236) largest angle weak lensing measurement of an individual cluster to date. 
Our analysis is performed using the SDSS which is the only imaging survey that covers a large 
enough area to measure the shear due to Coma. We find the shear can be fit well to an NFW 
profile and have compared our weak lensing derived mass estimate to dynamical methods 
which have been used to probe the mass of Coma out to large radius. In part icular we have 
comp ared our weak lensing mass to the mass derived using the virial theorem ([The fc White 
19861 ). X-ray data (lHugheslll989l ). and the cluster infall region (IGeller et al.lll999l ). We find 
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the virial mass of Coma from weak lensing is consistent with the mass derived using these 
other techniques, though the error from weak lensing is larger. 
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Table 1. Virial Mass Estimates of Coma 



Virial Radius 
(h _1 Mpc) 



Mass 
;iO 15 h" 1 M ) 



Reference 



1 qq+0.21 
±.yy_ 22 


i OO+0.65 
-■-■00-0.56 


this work a 


1.5 


0.8 


l a 


2.5 


0.93 ±0.12 


2 


2.7 


0.95 ±0.15 


3 



a Determined from the NFW r2oo virial 
radius. 



References. 



(1) 


Geller et al. 


1999 


The k Whitel 


1986. 
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Fig. 1. — r band magnitude distribution of source galaxies used in our analysis. Objects with 
extinction corrected model magnitudes between 18 < r < 21 are used. The total number of 
source galaxies used to measure the shear due to Coma is ~ 270, 000. 
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Fig. 2. — Tangential shear centered on the Coma Cluster in the SDSS. The measured shear 
is shown as solid squares along with la error bars. The 45° component is shown as open 
triangles and is consistent with zero. The solid line represents our best fit NFW model. 
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Fig. 3. — Average tangential shear measured from "blank" regions of the SDSS North not 
associated with Coma. The shear signal (solid squares) and 45° component (open triangles) 
should be zero in these regions. Results are the inverse variance weighted average of 6 
independent ~ 20° x 20° patches. The size of the annuli here match those used to evaluate 
the shear due to Coma. 



